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For ITER-relevant runaway electron studies, such as suppression, mitigation, termination 
and/or control of runaway beam, obtaining the runaway electrons after the disruption is 
important. In this paper we report on the first achieved discharges with post-disruptive 
runaway electron beam, entitled “runaway plateau”, in the COMPASS tokamak. The 
runaway plateau is produced by massive gas injection of argon. Almost all of the disrup¬ 
tions with runaway electron plateaus occurred during the plasma current ramp-up phase. 
Comparison between the Ar injection discharges with and without plateau has been done 
for various parameters. Parametrisation of the discharges shows that COMPASS disrup¬ 
tions fulfill the range of parameters important for the runaway plateau occurrence. These 
parameters include electron density, electric field, disruption speed, effective safety fac¬ 
tor, maximum current quench electric field. In addition to these typical parameters, the 
plasma current value just before the massive gas injection surprisingly proved to be 
important. 


1. Introduction 

As the tokamak concept developed in the last 50 years and advanced towards the 
ITER design, numerous challenges occurred and many were solved. One of the remaining 
tasks is control or mitigation of Runaway Electrons (RE) in ITER after the disruption. 
Estimations from codes predict RE with several tens of MeV to carry up to 70% of pre- 


disruptive plasma current (Hender et al. 2007 p. S178). As deposition of runaway electron 
beam can be highly localised, it could severely damage plasma facing components and 
blanket modules of ITER. 

The electron is said to run away, when the collisional drag force acting on it becomes 
smaller than the accelerating force coming from the toroidal electric field E tor . There 
are three main mechanisms for the runaway generation: a) Dreicer (primary) mechanism 
(Dreicer 1959| 1960); b) hot-tail mechanism (Smith & Verwichte 2008); c) avalanche 
(secondary) mechanism (Rosenbluth & Putvinski |1997 ). However, there is a theoretical 
limit for the electrical field, so called critical field E cr u , under which RE cannot be 
produced by these mechanisms (Connor~fe Hastie|~1975). The toroidal electric field E tor in 
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ITER during the stable discharge will be under the E cr i t threshold, making the controlled 
ITER plasma void of the RE. On the other hand, if disruption occurs, the electron 
temperature T e would drop during Thermal Quench (TQ), and thus plasma electric 
resistivity 77 would increase. E tori being proportional to r) j, will rise dramatically during 
the Current Quench (CQ), because the current density j drops much slower than the 
electric resistivity increase due to the vessel electromagnetic field penetration time. This 
increase of the field will first induce runaway seeds that will then be multiplied enormously 
by the avalanche effect. In the ITER disruption scenarios, the avalanche multiplication 
factor could be as large as 10 22 (Hen der et al. | |2007, table 5), forming an electron beam 
that could threaten ITER’s first wall structure. Following the above outline, ITER should 
be equipped with a proper suppression and/or mitigation technique dedicated to the RE 
control. Thus, achieving post-disruptive RE beam is one of the first significant steps for 
COMPASS towards the ITER-relevant runaway suppression/mitigation studies. 


The COMPASS tokamak (Panek et a/.|2 006) is a small-size experimental fusion device 
with major radius Rq = 0.56 cm and minor radius a = 0.23 cm. Toroidal magnetic field 
B tor is in 0.9 — 1.25 T range and plasma current I p can reach up to 330 kA. Electron 
densities are flexible and are typically of order of magnitude of 10 19 — 10 2 O m~ 3 . Plasma 
shaping varies from circular and elliptical to single-null D-shaped ITER-like plasmas. The 
typical pulse length is 0.4 s, although the low current circular discharge with RE can last 
almost 1 s. Furthermore, flexibility of various plasma parameters (e.g. shaping, densities, 
plasma current, etc.) combined with significant, but still safe, runaway population make 
COMPASS suitable for runaway models validation and scaling towards ITER. 


In contrary to large tokamaks (e.g. JET, ITER), where most of the RE are produced 
during the disruption ( Martin et a/.||1995| |Yoshino et aL|[l999 Gill et a7T||2000 1 , in small 
and medium size tokamaks RE are created either during the current ramp-up or the 
flat-top phase (Esposito et al. 2003 Paprok et al. 2013) when n e is low and/or E tor is 
high enough. Additionally, the present COMPASS maximum value for B tor is 1.25 T, 
while various observations noted that getting the post-disruptive RE spontaneously is 
not possible if B tor is under « 2 T (Martin et al. 1995[ Yoshino et al. 1999 Gill et al. 


2002). The R tor -limit is the most probable reason for the lack of post-disruptive runaway 
observations in COMPASS. Therefore, size of the COMPASS and its maximum B tor could 
discard this facility from the ITER-relevant runaway suppression/mitigation research. 
Nevertheless, some of optimism can be found in experiments in which the post-disruptive 


RE were achieved with high-Z Massive Gas Injection (MGI) (Yoshino et al. 

Hollmann et al. 2013) or high-Z pellet injection (|Yoshino et al. 1999| 


et al. 


et al. 


2002 


1999 Gill 


Hollmann 


2013), since some of these experiments had B tor lower than 2 T. Moreover, a detailed 


study of Rt or -limit as a function of amount of Ar injected was performed recently in 
JET (Reu x et qZ.| [2014), where post-disruptive RE were observed even for B tor = 1.2T. 
Therefore, Ar injection was used to trigger the first post-disruptive RE in COMPASS. 


The paper is organised as follows: in section [2j the experimental setup used for the 
experiments and demonstration of runaway plateau observation is presented. In section 
[3j general runaway parameters are reported, followed by the injection and disruption 
details. The section is finalised with the discharge analyses of the parameters important 
for the plateau occurrence. In section [4j the results presented in section [3] are discussed. 
Finally, in section [5] conclusions and future perspectives are addressed. 
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Figure 1: Principal diagnostics used for the RE plateau studies. 


2. Introduction to Experiment 

2.1. Experimental Setup 

In all discharges described in this paper the plasmas were circular - limited by the car¬ 
bon High Field Side (HFS) wall - with additional carbon Low Field Side (LFS) limiter 
for the inner wall vessel protection (see Fig. [lj. Typical magnetic field B tor was 1.15 T 
and plasma currents I p at the moment of the gas injection varied from 40 to 140 kA. 
The electron density n e was relatively low (0.8 — 2.2 x 10 19 m~ 3 ) for maximising the 
runaway generation. Schematics of the experimental setup used in the experiments pre¬ 
sented in this article is shown in Fig. |T] In this article, measurements of runaway losses 
will be presented from PhotoNeutron (PN) detector located nearby the north wall and 
Nal(Tl) scintillator for Hard X-Ray (HXR) detection located at the south-east part of 
the tokamak hall. Both detectors are approximately 5 m from the vessel. Photoneutrons 
with energy of few MeV are observed with the ZnS(Ag) neutron detector embedded in a 
plastic matrix. Beside the neutrons, the PN detector is suspected to be sensitive to the 
strong fluxes of HXR, although the detector is shielded by 10 cm of Pb. HXR are mea¬ 
sured with unshielded Nal(Tl) scintillation detector, where the signal is amplified with a 
photomultiplier tube and the energy range is approximatelly from 100 keV to few MeV. 
Furthermore, the low energy photon radiation measurements will be presented from H a 
detector and the bolometry. H a detector is located radially at the east part of the toka¬ 
mak vessel. AXUV photodiodes, located at the north-west part of the tokamak vessel, 
with photon energy response from 7 eV to 10 keV are used for bolometric measurements. 

MGI of argon was achieved using a solenoid valve, located on the east side of the 
tokamak. The solenoid gas valve is connected to the vessel through two stainless steel 
tubes: the first one is 20 cm long and has an inner diameter of 4 mm, while the second one 
is 40 cm long and has an inner diameter of 6 mm. This non-negligible tube length implies 
a delay between the time of valve opening and the time at which the argon puff starts to 
interact with the plasma, i.e. roughly the time at which the gas enters the vacuum vessel. 
The delay is estimated to be approximately 1 ms taking into account a mean velocity of 
approximately 400 m/s for argon gas in vacuum at 300 K. 
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The Ar flow rate dN/dtt through the injection system was evaluated experimentally as 
a function of the back pressure Pback an d with linear dependence as follows: 


— = (9.5511 pbacfe - 1.00 83) 10 2 °, (2.1) 

where pback is in bars and dN/dt is in particles/s. The pressure pback used for the 
plateau discharges were 2.4 and 1.2 bar, corresponding to particle flow rates of 2 x 10 21 
and 10 21 particles/s, respectively. The valve is roughly estimated to be open 2 ms, better 
knowlegde of gas valve performance will be available soon by the installation of a fast 
opening and more reliable valve. 

Since there is a pipe between the valve and the tokamak vessel, the duration time 
for Ar to enter the vessel is larger then the opening time of the valve and this Ar puff 
duration in the tokamak vessel will be estimated here. For the two aforementioned back 
pressures, 2.4bar and 1.2 bar, the manufacturer gives a flux through the solenoid valve 
at standard conditions of 100 and 50Pam 3 /s corresponding to flow rates of 2.4 x 10 22 
and 1.2 x 10 22 particles/s at 300 K, respectively. Therefore, assuming a constant flow rate 
of the solenoid gas valve with the increase of pressure in the stainless steel pipe and 
neglecting the flow rate through the injection system in the tokamak, one can calculate 
the number of particles that fill the pipe and that will be puffed in the tokamak later on. 
Remembering that the valve stays open for about 2 ms, one can find that there will be 
about 5 x 10 19 and 2.5 x 10 19 particles for 2.4 bar and 1.2 bar, respectively. Notice that 
these numbers are much smaller than the total number of particles that can be stored 
in the pipes at 2.4bar and 1.2bar (8 x 10 22 and 4 x 10 22 respectively), justifying the 
assumption of constant flow rate through the solenoid valve. Now, knowing the flow rate 
through the injection system and the number of particles in the pipes, one can give an 
estimation of what the Ar puff duration in the tokamak vessel is: 25 ms for 2.4 bar and 
12.5 ms for 1.2 bar. The runaway plateau created in this manner lasted from 2.5 to 10 ms. 


2.2. Plateau Observation 

An example of a typical COMPASS discharge with MCI generated runaway plateau is 


shown in Fig. 2a together with slow I p decay for comparison in Fig. 2b Fig. 2a shows the 


plateau discharge #8585, when Ar puff starts to cool down the plasma, I p starts to drop 
and plasma radiation increases. After approximately 2 ms (this delay will be justified in 
the next section), TQ occurs and almost all plasma energy is radiated. At the same time 
the HXR measurement shows relatively low peaks in half saturated state and PN signal 
is rather low, meaning that high energy RE created during the discharge initial phase 
are still confined. Then, during the CQ, E tor is increased and boosts runaway production 
creating and amplifying the runaway beam. After the CQ the runaway beam carries non¬ 
zero current called Ire, lasting for few milliseconds as one can see from the top graph 
in Fig. [2a[ Finally, the RE beam terminates with the loss seen in HXR and PN signals, 
while there are almost no H a radiation and radiated power P ra d from plasma proving 
the runaway plateau existence. On the other hand, the COMPASS discharge #8668 


see Fig. 2b I displays an example of the slow radiative decay with MGI on COMPASS 


resembling to I p ramp-down, for which no TQ and CQ (a typical sign of fast disruption) 
are observed. In this discharge plasma radiates on a long time scale (« 20 ms). Although 
HXR and PN signals show presence of released RE, we shall not consider this as the 
runaway plateau, because the I p current is mainly driven by the thermalised plasma 
and not by the runaway beam, as one can see from strong H a emission. Notice that 
the difference in H a and P ra d measurements makes the distinction between the runaway 
plateau and the slow radiative /„ decay. The former one has a relatively low radiation 
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(a) #8585 (b) #8668 


valve Ar in valve Ar in 



Figure 2: Time evolution of the COMPASS discharge #8585 as an example for runaway 
beam (a) and #8668 as an example for slow I p decay (b) both initiated by the MGI. 
Plasma current I p , electron density n e and loop voltage Vi oop are plotted on the top, in 
the middle HXR and PN signals are showing RE losses on the wall, while the H a and 
Prad, measurements are showing the radiation losses from plasma on the bottom (N.B. 
The y-axis are different for the bottom plot of the two discharges). 


level after the disruption, indicating that there is only cold plasma beside the runaway 
population. 

Supplementary to the previous description of the RE plateau, the observation of RE 
beam with visible camera is displayed in Fig(3] for discharge #8585. The creation and 
localisation of the beam are well visible. 


3. Results 

Out of 137 discharges performed during the COMPASS RE campaign, Ar puff was 
used in 39 discharges where only 5 discharges ended in spontaneous disruption. Out of 
the 39 discharges, 14 had the RE plateau after the Ar puff, while 9 resulted in slow 
radiative I p decay, similar to a ramp-down. The remaining 11 discharges ended in a 
typical COMPASS disruption, i.e. without any RE. 

Based on these observations, all discharges with the Ar puff can be classified as: 

(a) STRONG (RE plateau) - I RE > 5kA 

(b) WEAK (RE plateau) - I RE < 5kA 

(c) SLOW (radiative current decay) - plasma current slowly decreases in the similar 
manner as a ramp-down phase 

(d) ZERO (RE plateau) - “typical” disruption for COMPASS with no RE remaining 
or generated after the disruption 

An example of each class is shown in Fig. [4j where Fig. 4b is a zoom of Fig. [4a] to 
emphasize the difference between weak and zero plateau measurements. Although these 
two cases could seem identical at the first sight, the PN signal confirms release of the RE 
after the disruption in the weak case (b) and their loss during the disruption in the zero 
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Figure 3: Visual observation of the RE beam with visible camera for discharge #8585: 
(a) before Ar reaches vessel, (b) formation of the RE beam on HFS, (c) RE beam and 
(d) RE beam drifts towards LFS. 


(a) All four classes. (b) Zoomed area from (a). 



990 1000 1010 1020 995 997 999 1001 1003 1005 

time [ms] time [ms] 

Figure 4: Classification examples: (a) #8672 for strong plateau, #8673 for weak plateau, 
#8668 for slow plasma current decay and #8677 as an example of disruption without 
RE surviving or produced; (b) Zoomed region from part (a) for better observation of the 
difference between weak and zero plateau, as well as photoneutron signal for comparison 
of weak and zeros cases. 
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Figure 5: Time evolution of plasma vertical Z and horizontal R — Rq positions for the 
COMPASS discharge #8585, associated with the plasma current I p barycentre. Positive 
values of R — Ro mean that the plasma I p barycenter is closer to the LFS and positive 
values of Z mean that the plasma I p barycentre is closer to the top of the vessel. 


case (d). This classification is very important, as it will be used from now on throughout 
the paper. We shall now present the main results of the RE COMPASS campaign. 

First, we have estimated the typical maximum runaway energy for the analyzed dis¬ 
charges to be 10 —15 MeV by taking into account the electron acceleration due to electric 


field with the synchrotron radiation losses only, as suggested in Martin-Soils et al. (2010); 


Yu et al. (2013). 


Second, since the RE before or during the disruption are more likely to be produced in 


the hottest center of the plasma (Gill et al. 2002), the remaining and newly produced post- 


disruptive RE may have more peaked radial current profile than the pre-disruption I p 
profile. The peaking represents localisation of the plasma current I p around the magnetic 
axis and can be expressed through the internal inductance U, which is calculated by the 
EFIT reconstruction ( Havlfcek fe Hronova|2010 ) at COMPASS. It was observed that the 
li value increases by only 5-45% in comparison to measurements on JET (Loarte et al. 


2011). In addition to the modest k rise, the normalised plasma pressure f3 n rises above 
1.5 for the same discharges and thus confirms that the overestimated /3 n as seen by EFIT 


(Vlaini6 ~ef gZ.||2014 ) is caused by the presence of RE. 

The inward motion (towards the HFS wall, negative R — Rq values in Fig. [5]) of the 
post-disruptive plasma, followed by its return towards the vessel center is in agreement 
with the TFTR (Fredrickson et al. 2015) and Tore Supra ( Saint-Laurent et aL||2009 ) 
observations. However, TFTR and Tore Supra feedback systems were able to stabilise 
the runaway beam, while presently in COMPASS the beam continues to shift outwards 
until its termination, as shown in Fig. [5] Note that the outward shift is also visible in 
Fig. i The vertical plasma position for the majority of cases is rather stable (an example 
being given in Fig. [5j , only in a few discharges some downward shifts were noticed. 


3.1. Disruption Generated by Argon 

As already mentioned, a solenoid valve was used to inject Ar gas into the plasma. Even 
though two different pressures were used (2.4 and 1.2bar), no particular differences in 
runaway beam parameters were identified. The reason could be that the pressure was 
only varied by a factor of 2. 

In devices larger than COMPASS, high-Z gas injection is used to trigger fast CQ in 
order to improve runaway generation (Reux et al. J2014]). The plasma current quench rate 
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Figure 6: Ioh as a function of the plasma current before the gas puff Idisr- Negative 
values of Ioh correspond to the current ramp-up phase, while positive values represent 
the current flat-top and ramp-down phases. 
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is the quantifying parameter for the CQ speed. The calculation of J 7 values for disruptions 
with and without Ar puff was performed. No particular differences were observed between 
the discharges, as the majority of the J 7 values are between 500 —2500 s -1 in both cases. 
The values implicate that the whole pre-disruptive I p is lost in 0.4 — 2 ms, which is 
the order of magnitude of the electromagnetic field penetration time of the COMPASS 
vacuum vessel (~ 0.5ms). 

3.2. Parametrisation of Runaway Plateau 

The ohmic heating (OH) central solenoid current Iqh - called MFPS in (Havlfcek & 


Hronova|2008 ) - will be used to indicate on the appearence time of the runaway plateau. 
For the RE discharges analyzed in the article Ioh is negative during the current ramp-up 
phase followed by Ioh at zero value for few milliseconds during the transition towards 
the current flat-top phase. For the rest of the discharge, i.e. current flat-top and ramp- 


down, it becomes positive and controlled by the feedback system (Janky et al. 2014). 


In Fig. 6} Ioh 2 ms before TQ is plotted versus I p also taken 2 ms before the TQ and 
denoted as Idisr- The reason why exactly 2 ms are taken will be seen later in this section, 
but it can be explained as the time before Ar starts to cool down the plasma, displayed 
in Fig. |4a| Also, later in the article, the measured parameters denoted with the index 
disr (e.g. Edi sr and ndisr) are taken at the same time. 

Fig. [6] shows that only one weak plateau out of 14 plateau discharges did not appear 
during the ramp-up phase, but during the ramp-down phase. Hence, RE plateaus are 
more likely to be produced during the current ramp-up phase than during the flat-top. 
The ramp-down case requires further investigation in future experiments, as only one 
such discharge was observed. 


Yosliino et al. (1999) did the first detailed parametrisation of disruptions with runaway 


occurence in JT-60U tokamak. According to his article, study of / 7 versus g e // is impor¬ 
tant for the plateau occurrence, where the effective edge safety factor q e ff for circular 
plasma is defined as: 


Veff = 


5 a 2 B t , 
Rio 


1+ (f? 


1 + 


(/3 p + U/2Y 


(3.2) 
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Figure 7: Plasma current quench rate / 7 as function of q e ff- Vertical red line corresponds 
to q e ff = 3.5 and horizontal blue line corresponds to / 7 = 500 s -1 . 
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Figure 8: Normalised electric field E disr /E crit (a) and plasma current just before the 
MGI puff Idisr (b) as function of the electron density n d i sr . The vertical blue line matches 
n d isr = 1-4 x 10 19 m -3 . The horizontal red line in (a) corresponds to E d i sr /E cr i t = 250, 
while black line in (b) is for I d i sr = 120 kA. The oblique green line in (b) represents the 
limit between the strong and the weak case. 
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The internal inductance and the poloidal beta /3„ are taken from the EFIT reconstruc¬ 


tion at the closest moment from the disruption. J 7 is already defined in Eq. 3.1 Fig. [7] 


shows / 7 versus g e // for the case of COMPASS. For all plateaus except the slow ones, J 7 
is between 500 and 1800 s -1 and q e ff is between 2.5 and 8. It is interesting to observe 
how majority of the zero disruptions are under q e ff = 3.5. Obviously, slow disruptions 
have significantly slower current decay than the rest of the discharges, their / 7 values are 
under 100 s -1 . 

According to the theory the production of the RE is more intense for lower densities. 
Thus, E disr normalised to E crit and I d i sr are plotted as a function of the line averaged 
density n disr measured by the interferometer in Fig. [8] Approximately, the critical value 
of electron density for obtaining the runaway plateau seems to be 1.4 x 10 19 m -3 . The 
ratio E d i sr /E cr u represents the relative strength of VJ oop . Critical value of the E d i sr /E cr n 
ratio on Fig.[8a]is around 250 for the analyzed discharges. Figure [8b| shows that the strong 
plateaus are created for lower I d i sr value than weak plateaus, taking the same n d i sr value. 
In addition, no strong plateau is observed above I d i sr = 120 kA, while half of the weak 
ones have I d i sr above 120 kA. 

Another parameter of interest is the current carried by the RE beam Ire- The depen¬ 
dence of Ire on the I disr is shown in Fig.[9aJ where only the ramp-up Ar MGI discharges 
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(a) (b) 
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Figure 9: Runaway beam current Ire as function of the pre-disruptive current Idisr (a) 
and the maximum loop voltage during the current quench Eqq (b). The colors are kept 
the same as in previous figures, while different marker type corresponds to different times 
of the Ar injection: 975ms (squares), 985ms (triangles), 995ms (diamonds). Symbols 
stand for a mean value of the Ire , and error bars stand for maximum and minimum 
values of the Ire- The vertical lines in (a) match 70kA (black), 100kA (blue) and 
120kA (red). 


are presented. The discharges in Fig. [9] are grouped by the time of Ar puff. Typical weak 
plateau Ire is between 0.5 — 3.5 kA, while Ire for strong plateaus decreases with Idisr 
and time of the puff. Furthermore, for the Ar injections performed at 985 ms and 995 ms 
the upper limit of the Idisr is indicated, namely 100 kA and 120 kA respectively. For val¬ 
ues lower than these Idisr values strong plateau seems to be produced, while above either 
weak plateau or no plateau occurred. In contrast to this, for the Ar injection at 975 ms 
the lower limit of Idisr is observed for about 70 kA, under which no plateau was detected. 
Anyway, more statistics are required. The dependence of Ire on the maximum electric 
field Eqq during the current quench (Fig. |9b[ ) has similar behavior like in Fig. 
could expect from the self-inductance effect between Idisr and Eqq- 


9a 


as one 


4. Discussion 

Even though the number of discharges devoted to runaway plateau studies was limited 
on COMPASS in the dedicated RE campaign, it was still possible to do comparative 
analyses. The results presented in the previous section are discussed in following order: 

• general characteristics on RE are outlined 

• observed differences between discharges with and without the runaway plateau are 
reported 

• issues on obtaining the RE plateau with Ar puff are discussed 

• achieving strong plateau is commented 

In Loarte et al. (2011), an increase of li by factor of 2 to 3 has been reported for the 


RE plateau, which is significantly larger then the modest rise observed in COMPASS 
(0.05-0.45). This modest k rise could mean that RE seeds generated during the ramp- 
up phase are occuring in larger relative area than for the case of JET. Similarly to Tore 
Supra ( Saint-Laurent et al.|2009 ), JET ( Loarte et al.|20ll ) and TFTR (Fredrickson et al. 
2015), the inward motion of the RE beam is observed for the beginning of the plateau 


phase at COMPASS. 

Almost all (13 out of 14) generated RE plateaus were achieved for Ar puff in the ramp- 


up current phase, as for Tore Supra (Saint-Laurent et al. 2011). Regarding disruptions 
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itself, the CQ speed is one order of magnitude larger than for the case of JT-60U (Yoshino 


et aL||l999 ), where I 7 > 100 — 200 s 1 was reported as the plateau formation condition. 


Concluding that COMPASS has fast enough disruption for the plateau formation (see 
Fig-0. but other factors - e.g. B tor , Vi oop , avalanching ( |Rosenbluth fe Putvinski|[l997| 
- are not fulfilled, explaining why the Ar MGI is necessary for COMPASS to obtain the 
runaway plateau. 

Next, q e ff and its relation with / 7 is one more important plasma characteristics for the 
plateau creation. In the case of JT-60U ( |Yosh in o et aZ.|1 999), on top of the / 7 condition, 
q e ff has to be over 2.5. From Fig. [7] it is obvious that disruptions analyzed here are 
deep in the reported parameter region. However, there is an indication how the plateau 
condition for / 7 and q e ff in COMPASS could be different from those observed in JT-60U. 
Anyhow, this is still to be investigated by enhancing the statistics. 

As observed from Fig. |8a| the limiting ndi sr for plateau to appear is around 1.4 x 
10 19 m -3 , which corresponds to E cr i t = 0.0122 V/m. The Dreicer mechanism is the most 
probable source of the post-disruptive production of RE at COMPASS, because the 


avalanching is expected to be important for the tokamaks with Ip > 1 MA (Rosenbluth 


& Putvinski 1997). However, from Fig. 8a it is apparent how the strong plateaus are 


obtained for low £J isr / E cr it values compare to the weak plateaus. This observation, at 
the first sight contra-intuitive, could be explained with the appearance of the avalanching 
effect, as avalanching is dominant runaway generation mechanism for lower Edi sr /E cr n 
assuming that the electron temperature profile remains unchanged (see Nilsson et al. 
2015 Fig. 10). Anyhow, this possibilty is still to be investigated. The Dreicer field Ep > 
is currently difficult to determine as the Ar injections were often too early, so that no 
Thomson scattering data were collected yet. For the cases plotted in Fig. [8bJ it seems 
that lower densities are necessary in order to achieve plateau for Idisr above 120 kA, 
making Idisr important parameter for the plateau production. 

In the COMPASS case, it seems that the inverse dependence is recognised for the 


strong plateaus (Fig. 9aI. The dependence of Ire on Idi sr from Fig. 9b looks almost the 
same as the one from Fig. [9a] as one would expect. This observation comes from the fact 
that the amplitude of induced E tor during the CQ is directly proportional to the I p before 
the disruption. Lower and upper boundary signs of plasma current for strong plateaus 


from Fig. [9b] are not unique, these boundaries have been observed in JET by Gill et al. 
(2002). In this article, the lower limit is assigned to low E tor , while the upper limit is 


possibly connected to magnetic fluctuations. For COMPASS more discharges would be 
required to improve the statistics and find the two limits. 


5. Conclusion and Future Work 

Before this dedicated campaign, runaway plateau was never observed in COMPASS. 
As a matter of fact, there was scepticism concerning the possible plateau occurrence for 
any plasma condition, due to the COMPASS tokamak size, low B tor and the low plasma 
currents leading to relatively low electric field E tor during the disruption. Nevertheless, 
in this paper a clear demonstration of obtaining runaway plateau by MGI is reported. 
The RE plateau currents varied between 0.5 to 40 kA, with duration from 2.5 to 10 ms. 

Argon injection disrupted discharges in COMPASS have been investigated in order to 
clarify the necessary conditions for runaway plateau production. It was found that the 
easiest way to produce the RE plateau was to inject Ar during the ramp-up of the plasma 
current. Furthermore, the typical COMPASS disruptions without RE can satisfy various 
parameters important for the runaway plateau creation (e.g. n e , V) oop , / 7 , q e ff) without 
Ar injection, thus high-Z MGI is probably required only for activating thermal quench 
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to enhance runaway population. Unusually, for the discharges considered in the paper, 
it seems that the plateau generation also depends on plasma current during the Ar puff 
injection. Even though, the CQ after MGI induced disruption has a very short time, it is 
possible that avalanche mechanism is present in COMPASS during the runaway plateau 
formation. 

More experiments need to be done in order to draw final conclusions on the definite 
conditions for the runaway plateau generation in COMPASS tokamak. From present 
knowledge we can conclude that some observations correspond to reports from larger 
tokamaks, although the amplitude is sometimes different. Indeed, this difference in the 
magnitudes could be important for scaling towards ITER. 

The experiment presented here confirm that COMPASS is a tokamak suitable for 
various ITER-relevant runaway studies, such as: 

(a) studies of runaway plateau termination - energy balances and timescales (Loarte 
| et al . | |2011| |Martfn-Soh's et aI7||2014[ ) 

( b) improvements of the runaway beam mitigation 

(c) testing the runaway control system 

( d ) benchmarking of the runaway models 

Nonetheless, the scenario for inducing the runaway plateau is necessary before further 
ITER-relevant studies are performed. Presently, LUKE (Dec ker fe Peysson[ [2004) code is 
being used in collaboration with CEA for a better understanding of the physics behind 
the measurements. 
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